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Abstract 

The constituents of cold dense matter are still far from being understood. 
However, neutron star observations such as the recently observed pulsar PSR 
J1614-2230 with a mass of 1.97 ±0.04 M help to considerably constrain the 
hadronic equation of state (EoS). We systematically investigate the influence 
of the hyperon potentials on the stiffness of the EoS. We find that they have 
but little influence on the maximum mass compared to the inclusion of an 
additional vector-meson mediating repulsive interaction amongst hyperons. 
The new mass limit can only be reached with this additional meson regardless 
of the hyperon potentials. Further, we investigate the impact of the nuclear 
compression modulus and the effective mass of the nucleon at saturation 
density on the high density regime of the EoS. We show that the maximum 
mass of purely nucleonic stars is very sensitive to the effective nucleon mass 
but only very little to the compression modulus. 
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1. Introduction 

The equation of state (EoS) in the vicinity of saturation density is pretty 
well understood. But beyond saturation, the theories of dense matter present 
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uncertainties. Neutron stars provide a fantastic astrophysical environment 
for testing theories of cold and dense matter. In the core of neutron stars, 
densities could reach values of several times 10 15 g cm -3 . At such high den- 
sities, the Fermi energies of the constituent particles could exceed the rest 
masses of heavier particles, and hence favour the appearance of such particles 
in the core. Further, as the timescales associated with neutron stars are much 
greater than those associated with weak interactions, violation of strangeness 
conservation due to weak reactions in the core would result in the appearance 
of strangeness-containing particles such as hyperons. By producing new de- 
grees of freedom, the appearance of strange particles is expected to result in a 
softer EoS of dense matter in the neutron star interior. The highest neutron 
star mass that can be supported depends crucially on the EoS. Observations 
of pulsars in neutron star binaries provide a precise measurement of neutron 
star masses through general relativistic effects. The best determined pulsar 
mass (1.4414 ± O.OOO2M ) is the Hulse- Taylor pulsar, and masses of most 
other pulsars are found to be clustered around this canonical value. Recently, 
several neutron stars with larger masses have been discovered. Radio timing 
observations of three Post Keplerian parameters led to the most precise mea- 
surement of the mass of a millisecond pulsar of 1.667 ±0.021M Q [lj. Shapiro 
delay measurements from radio timing observations of the binary millisec- 
ond pulsar PSR J1614-2230 indicated a mass of 1.97±O.O4M of the neutron 
star [2]. This is the largest rather precisely observed pulsar mass so far and 
thus poses the tightest reliable lower bound on the maximum mass of neu- 
tron stars. Any theory of ultradense matter requires that the EoS produce 
a maximum mass at least as high as this measured value, i.e., models with 
M mm (theo) < M max (ohs) would be ruled out. 

According to existing models of dense matter, the presence of hyper- 
ons leads to a considerable softening of the EoS, resulting in a reduction of 
the maximum mass of the neutron star. With hyperons, including only the 
hyperon-nucleon interaction, Brueckner-Hartree-Fock (BHF) calculations ob- 
tain maximum masses of the order 1.47M j3|]. The inclusion of the hyperon- 
hyperon interaction leads to a further softening of the EoS and reduces the 
obtained masses to 1.34M Q Q. On inclusion of three-body interactions, the 
maximum mass achieved in the BHF framework is 1.26M By employ- 

ing a recently constructed hyperon-nucleon potential, the maximum masses 
of neutron stars with hyperons is computed to be well below 1.4M C 



other approach is to adopt a relativistic mean field (RMF) Model [8l l9l. llO 



which we employ for the following investigation. Several attempts have been 
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made using the RMF model to explain neutron star masses higher than 2M Q , 
by artificially increasing the hyperon vector coupling away from their SU(6) 



values [ll|, [12|, [13J. Dexheimer et al. |14| obtained large neutron star masses 



using a chirally motivated model including fourth-order self-interaction terms 



of the vector mesons u, p and <f>. Recently, Bednarek et al. 15j achieved a 
stiffening of the EoS by invoking quartic vector-meson terms proportional to 
o; 4 ,p 4 ,0 4 and cross terms like <fr 2 uj 2 ,<j) 2 p 2 and u 2 p 2 . When the stiffening on 
inclusion of hyperons is not sufficient, a transition to the quark phase can be 
considered to produce large maximum masses [l6j]. Bonanno and Sedrakian 



17l | also succeeded in obtaining a large neutron star mass with a hyperon 



and quark core using a stiff EoS and vector repulsion among quarks. 

As the parameters in the RMF model are fitted to the saturation proper- 
ties of infinite nuclear matter, extrapolation to higher densities and asymme- 
try involve uncertainties. Three of these properties - the saturation density, 
the binding energy and the asymmetry energy are more precisely known than 
the remaining ones - the effective nucleon mass and the compression modu- 
lus of nuclear matter. The uncertainty in the dense matter EoS is basically 
related to the uncertainty in these two saturation properties. In this paper, 
we determine which of the two influences most the high density part of the 
EoS that determines the highest attained neutron star mass. The parame- 
ters associated with attractive interaction among hyperons and nucleons are 
fitted to the potential depths of hyperons in nuclear matter, known from 
hypernuclear experiments. We investigate in this paper how the uncertainty 
in hyperon potential depths influences the stiffness of the EoS and hence the 
maximum mass. We want to address the question: under which conditions 
can hyperons be present in massive neutron stars? We achieve this through 
a controlled parameter study within the RMF Model. In our model, we suc- 
ceed in generating stiff EoS by taking small values of effective nucleon mass 
m*, in the range of those obtained from fits to nuclei rather than bulk nu- 
clear matter, and by including the strange vector meson 0. In a subsequent 



paper [18], we will focus on the assumption of the underlying symmetries 
that govern the repulsive interactions among hyperons and nucleons. 

This paper is organized in the following way: In Sec. 2, we describe the 
model to calculate the EoS. The parameters of the model are listed in Sec. 
3. The results of our calculations are discussed in Sec. 4, and the summary 
and conclusions are given in Sec. 5. 
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2. Theoretical Model 



One of the possible approaches to describe neutron star matter is to adopt 
a RMF model subject to chemical equilibrium and charge neutrality. For our 
investigation of nucleons and hyperons in neutron star matter we will choose 
the full standard ( J p = | + ) baryon octet as well as electrons and muons. In 
this model, baryon-baryon interaction is mediated by the exchange of scalar 
(<rV vector (cj) and isovector (p) mesons. The Lagrangian density is given 
by 

C = ^2^ B -m B + g aB (T - QuBl^ ~ g P Bl^B ■ P M ) ^B 

B 

+ i (d^a - m 2 a a 2 ) - U(a) + U(u) 

+ -m 2 u u^ - ^ Pfll/ ■ p» v + -m 2 pPfl ■ pP. (1) 

The isospin multiplets for baryons B are represented by the Dirac spinor b 
with vacuum baryon mass tub, isospin operator t#, and and p^ u are field 
strength tensors. To reproduce the saturation properties of nuclear matter, 
the scalar self-interaction term 

U(a) = V + lea 4 (2) 



is introduced [19|. We also included an additional self- interaction term 

U(u) = \d{uj^f (3) 



for the vector field as proposed by Bodmer [20] . Due to inclusion of this term, 
the vector field increases proportional to p 1 / 3 for high densities, where p is 
the baryon density, instead of the linear dependence in absence of this term. 
This results in a good agreement with Brueckner-Hartree-Fock calculations. 
We will denote the above model with the three exchange mesons as "model 
acup" . The hyperon-hyperon interaction is usually incorporated through the 
exchange of additional strange scalar (cr*) and strange vector (0) mesons, 
again the scalar meson being responsible for attractive and the vector meson 
for repulsive interactions respectively: 

C Y Y = Yl ^ B (9-*BV* ~ g+BlpV) ^B 
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~^r v + ImfaP. (4) 

Since it is our goal to obtain the stiffest possible EoS within the model, we 
will only make use of the <fi but omit the a* . This choice is also in accordance 
with AA-hypernuclear data, where now it is clear that the AA interaction is 



only weakly attractive (see Ref [21[ for a discussion). The model obtained 
with this additional Lagrangian we will call "model acopcp" . The calculation 
is performed using the mean field approximation (22|. The effective baryon 
mass is given by m* B = itlb — g a B<J while the chemical potential of baryon B 
is /ie = {kp B + m*B ) 1/2 + 9wBUo + 94>b<Po + h B g P BPm/^^ with k Fg denoting 
the corresponding Fermi momentum. Charge neutrality is described by the 
condition 

Q = ^ QBn B ~ n e - = , (5) 

B 

where ub is the number density of baryon B, qs is the electric charge and n e 
and n M are charge densities of electrons and muons respectively. We consider 
the hadronic matter to behave like an ideal fluid. The energy-momentum 
tensor for such a fluid yields the EoS, defined by the relationship between 
the total energy density 

e = X -mla 2 + ha 3 + + 

1 2 2,1 2/2,1 22 

+ 2 m ^ w o + + 2 m pP03 

+ ^ 2 -^l kFB (k 2 + mB 2 ) l/2 k 2 dk 

E ^ I (k 2 + m 2 Y /2 k 2 dk, (6) 



l=e~ ,ji 



and the pressure 



P = -\mlo 2 -\bo* -lea* + \d^ 

,1 2 2,1 2/2,1 22 

+ 2 m ^o + 2 m A + 2 m pPos 
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1 ^ 2J B + 1 f kFB k 4 dk 
3 \ 2vr 2 J (k 2 + mfy/ 2 

'3 J-^ vr 2 i (A; 2 + m 2 ) 1 /2 ' ( ^ 



3. Parameters of the model 

5.1. Nucleon-Meson coupling constants 

There are five input parameters for the RMF models: the nuclear satura- 
tion density p ( or equivalently, the number density n = po/m N ) as well as 
the binding energy per baryon number B/A, the effective mass of the nucleon 
m* N , the nuclear compression modulus K and the asymmetry coefficient a sym , 
all taken at saturation density. The parameters g^, g^N, 9 p n, b and c of 
"model aup" are determined from the saturation properties of nuclear mat- 



ter in an analytic way (see [23J for explicit formulae). The binding energy, 
the asymmetry coefficient and the saturation density are well determined. 
We set them to n = 0.153 /m" 3 , B/A = -16.3 MeV, a sym = 32.5 MeV (2 



The effective nucleon mass and the compression modulus are less well known. 
In principle, the parameters of the model can be fitted to properties of nuclei, 



which result in low effective masses [25j. However, here we want to explore 



the full parameter range of the model, as model parameters fitted to bulk nu- 
clear matter and usually adopted in the community have large effective mass 
parameters. Hence for this study, we use the parameter sets such as GM1 
(K=300 MeV, m* N /m N = 0.7) and GM3 (K=2A0 MeV, m* N /m N = 0.78) 
which are fitted to the properties of bulk nuclear matter and have rather 
large effective nucleon masses in the range m* N /mN ~ 0.7 — 0.8 j24j], and also 
the parameter sets fitted to the properties of nuclei, e.g. the whole NL pa- 
rameter family and TM1 which yield very low values of m* N /mN ~ 0.5 — 0.65 



26|,|27 



The value of compression modulus for symmetric nuclear matter can be 
extracted from the energies of the isoscalar giant monopole resonances (IS- 



GMR) in nuclei (see e.g. [28] and references therein). Systematic studies 
of ISGMR energies in various nuclei suggest a value of 231 ± 14 MeV for 
the incompressibility of symmetric nuclear matter (29I . 3o| . However, there 



are strong surface effects, as well as large uncertainties in the value of the 
compression modulus for higher densities and asymmetric nuclear matter. 
Another semiempirical parameter of nuclear matter crucial for determining 
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the stiffness of the EoS is the density dependence of symmetry energy, de- 
noted by L. We computed the values of the density dependence of symmetry 
energy L for the sets discussed above, and they were found to be fairly com- 
parable with the values existing in the literature 



31. 32 



These however 

differ from the values of L extracted from isospin diffusion data (L = 88 ± 25 
MeV) [33[ and isoscaling data {L ~ 65 MeV) We note here that when 
the effective mass m*/m is varied between 0.55 and 0.8, L changes from 108 
to 94 MeV (± 1 MeV depending on the compressibility). 

3.2. Hyperon-Meson coupling constants 

The hyperon coupling constants q^y, q P Y and g$y are determined by using 



SU(6) symmetry (the quark model) [35l. 136 



g P N — -g P T, 



9 P A 



2#0A — 2g0 S 



9uB 

0, 

2^2 



(8) 



The scalar meson-hyperon coupling constants g a y are adjusted to the poten 



tial depths Uy felt by a hyperon Y in a bath of nucleons at saturation [24 
following the relation 



(TV) 
Y 



9aYO- q + guY^Q 



eq 



(9) 



where the hypernuclear potential depths in nuclear matter are fixed in 
accordance with the available hypernuclear data. The best known hyperon 



potential is that of the A, having a value of about = —30 MeV [37J, [38, 



39] . In case of the Ss and the Ss the potential depths are not as firmly known 



as for the A. The following values U^ N) = +30 MeV, ui N) = -18 MeV are 



generally adopted from the hypernuclear experimental data (see I40j for a 
discussion). 



4. Results and Discussions 

4-1. Hyperon Potential Depths 

We begin our investigations with "model aup" , namely with a a—u model 
including non-linear scalar and vector self-interactions as well as p mesons 
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and the full baryon octet. In the previous section, we already mentioned that 
there are uncertainties in the values for the potential depths of hyperons 
in nuclear matter, obtained from hypernuclear experiments. We vary the 
potentials and systematically to study how they influence the 

stiffness of the hadronic EoS and the limiting neutron star mass. 

4-1.1. Varyin g 



We start with a fixed 5 potential depth of ui N) = -28 MeV [10|. The 
potential depth of E is varied in the range -40 MeV < U^ N) < +40 MeV. We 
plot in Fig. [IJthe EoS for different values of (lower branch) for the GM1 
parameter set and also for the same model with additional 4> meson (upper 
branch). For the sake of clarity only the cases = —40, —30, —20, 0, +40 
MeV are displayed. 

It can be inferred from the figure that for less deep potentials the 
EoS becomes stiffer for each model. However, this is only true for negative 
potentials, since the values U^ N) = MeV and U^ N) = +40 MeV basically 
result in the same EoS. For negative values of the potential the Ss are 

bound in nuclear matter and the deeper the potential is, the more attractive 
must be the effective mesonic interaction and thus the softer the EoS. For 
U^ N) > MeV, the Ss are no longer bound to nuclear matter and the ef- 
fective mesonic interaction becomes more and more repulsive for increasing 
U^K This would in principle stiffen the EoS. However, up to neutron star 
densities, i.e. about n B < (4 — 7)n , the Ss are not present in hadronic 
matter if the potential is repulsive and thus the EoS up to these densities 
becomes insensitive to the actual value of C/^ . The composition of hadronic 
matter for positive and negative values of £ potential depth has already been 
investigated in 

The repulsive potentials > give the highest maximum star masses 
for each parameter set (also for others which we did not plot, as e.g. TM1 
and NL3). In Fig. [2] we plot the parts of the mass-radius relations around 
the maximum mass star configurations for the EoS in Fig. [TJ The maximum 
mass for the lower branch of curves in Fig. [2] corresponding to "model croup 11 
is M max = 1.75M . We note, that the variation of [/^ in "model croup 11 for 
GM1 does not help to increase the maximum neutron star mass above the 
limiting value of 1.97±O.O4M . Investigating other parameter sets like TM1 
would not solve the problem since the maximum masses are comparable to 
those of GM1. The various NL-models would be good candidates, since they 



8 




Figure 1: EoS for different potential depths in MeV, indicated by the corresponding 
numbers in the key. The lower branch is obtained for "model cajp" for the GM1 parameter 
set while the upper branch is obtained including additionally the 4> meson. 
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Figure 2: Mass radius relations for neutron stars obtained with the EoS from Fig. Q] The 
variation of in "model aujp LL cannot account for the observed neutron star mass limit 
(lower branch), unless the </> meson is included in the model (upper branch). 
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are amongst the stiff est available EoS [10]. However, the NL- models have 
(like all models which are fitted to finite nuclei) very small effective nucleon 
masses already at saturation density (m*/m ~ 0.55 — 0.65). At around five 
times nuclear saturation density the a— fields are so large that the effective 
nucleon mass becomes zero or in principle even negative. This causes a well- 
known instability which makes it impossible to find a physical solution to 



the RMF equations [10[. However, the critical density for the onset of the 
instability is higher than the maximum central density in the neutron star 
interior, and hence this is not a relevant problem for the present investiga- 
tion. 

The "model aupcf>" gives us the upper bundle of EoS in Fig. [TJ In Fig. [2] 
the upper branch of the mass-radius relations corresponds to these EoS from 
"model cxa;p0". While the variation of £ potential only varies the maximum 
mass in the range Am w 0.03 M Q , the impact of <fi on the maximum mass is 
an order of magnitude larger: the maximum mass is increased by nearly 0.21 
M Q to M max =1.96 M & . We should also mention that for purely nucleonic 
matter EoS the obtained maximum masses are even higher (about 2.4M for 
GM1 and nearly 2.8M for NL3). 

4-1-2. Varying U~ N ^ 

We have seen that the potential depth has but a little influence 

on the maximum mass of a hadronic neutron star. As it does not play an 
important role in our investigations, we shall fix it to the value U^=+30 



MeV [lOj which means that the Ss are not present and the EoS is as stiff 
as possible. We then vary the S potential in the range —AOMeV < 

< +40 MeV. As before, we compute the EoS for the GM1 parameter set 
for "model aojp" and "model aojpcj)" for different H potentials. We plot in 
Fig. El the results for the values of from -40 MeV to +40 MeV in steps 
of 20 MeV. 

We see that for increasing H potential, the EoS stiffens. A difference to 
the influence of the S potential depth is found for positive values of L^i : 
here the EoS remains sensitive to the value of the potential due to the contri- 
butions from S~ at neutron star densities while for even higher densities both 
Ss contribute. In Fig. HJ we display how the composition of the neutron star 
core varies for various values of S potential in models u aup" and u aup4>" . 
It is clear from the figure that for positive S potential, A hyperons dominate 
particle fractions at high densities, while for negative the particle frac- 
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Figure 3: EoS for different potential depths l/i^ . The lower bundle of curves corresponds 
to "model aujp" and the upper to "model aujpcj)" . All EoS are obtained for GM1 parameter 
set. 
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Figure 4: Particle fractions for "model crwp' ' for potential depths f/4 = + 28 MeV ( a ) 
and -28 MeV (b) are displayed. Similarly, particle fractions for "model auipcf)" are shown 
for S potentials +28 MeV (c) and -28 MeV (d). The threshold for appearance of the S 
hyperons is pushed to higher densities with increasing 3 potential. 
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tions from all hyperons contribute. As S potential increases from negative to 
positive values, the threshold density for appearance of S~ and S° hyperons 
is pushed to higher densities, and this effect further increases on inclusion of 
repulsion from the (f> mesons. 

In Fig. [5] we plot the mass-radius relations of the neutron stars ob- 
tained from the various EoS in Fig. [3J Within "model aup" the maximum 
masses range from 1.69 M Q for ui N) = -40 MeV to 1.92 M© for ui N) =+A0 
MeV, i.e. over an interval of Am M . The maximum mass for the 

case L r | 7V) =+40 MeV including 4> is 2.04M Q . Although the value of the H 
potential depth helps to stiffen the EoS far better than the it is still not 
enough to be consistent with the mass of PSR J1614-2230 for "model aup" . 
However, the mass constraint can be fulfilled within "model aupcft" for the 
GM1 parameter set. 

4.1.3. The U^ N) - ujP Plane 

We combine the variations of and and plot in Fig. [6] lines 

of constant maximum mass in the C/g — U% plane within the range —40 
to +40 MeV for both potentials. As before, we compare the results from 
"model aup" with those of "model aupcfi". In the light of our previous 
discussions we can quickly explain the shape of the lines in the plot: On 
each side of the C/^ = £/| line the maximum masses are dominated by 
that potential which is smaller and are rather insensitive to the value of the 
larger potential. This means, that in the region where U^ N) > U^ N) the 
lines of constant maximum mass are parallel to the £7^ — axis. In the other 
part of the plot, where Z7jj < the behaviour is mirrored. Since the 

influence on the maximum masses is larger for smaller potentials, the mass 
lines (which are plotted for equidistant mass steps of O.O1M ) are densest 
towards the lower and the left parts of the plots but thin out with increasing 
potentials. We conclude this section by noting that for the whole range 
-40 MeV < U^\ui N) < +40 MeV the maximum masses of "model aup" 
for the GM1 parameter set are not compatible with the observed pulsar 
mass 1.97±O.O4M . Inclusion of the <fi meson helps to sufficiently stiffen 
the EoS and the whole range of -40MeV < U^ N) ,ljt N) <+40MeV is now 
allowed. However, in this case the new mass limit does not help to restrict 
the parameter range of the hyperon potentials. 
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Figure 5: Mass-radius relations obtained from the EoS of Fig. [3] As before, the variation of 
U~ N ^ in "model trwp" cannot explain the observed neutron star mass limit (lower branch), 
unless the <f> meson is included (upper branch). 
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Figure 6: Lines of constant maximum mass (key numbers in M©) in the — U~ N ^ 

plane. Upper plot for "model atup" in GM1 parameterization, lower one for the same 
model including also <f> mesons. 
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Figure 7: Maximum masses of neutron stars as functions of the effective nucleon mass 
for different compression moduli using "model crwp". The upper curves refer to purely 
nucleonic stars while for the lower curves hyperons are included. 
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4-2. Compression Modulus and Effective Mass 

In the previous section, we found that the variation of hyperon potentials 
were not sufficient to explain the observed large neutron star mass. We then 
proceed to find out how the maximum mass depends on the nucleon cou- 
pling constants. We vary the least well known of the saturation parameters, 
namely the effective nucleon mass m* N and the compression modulus K. We 
first consider purely nucleonic stars, and then repeat the exercise including 
hyperons in the "model aup" . In Fig. [7] we plot the maximum masses as a 
function of the effective nucleon mass at saturation m* N /mN for several values 
of K. We see that the maximum masses for the nucleonic case (upper bunch 
of curves) show a strong dependence on the effective nucleon mass: they drop 
from 2.88M at m* N /m N = 0.55 to 1.87-2.O1M at m* N /m N = 0.8. From 
Fig. Hit is clear that for moderate effective masses, m* N /m^ « 0.6 — 0.7, the 
compression modulus gains a little influence on the maximum mass for pure 
nucleonic stars, and causes a mass difference of less than O.O5M . Above 
this value of the effective mass, the compression modulus causes a larger 
splitting of the maximum masses until at m* N /mN = 0.8 the mass difference 
is about O.1M . For the hyperonic case (lower bunch of curves), the maxi- 
mum masses range from 2.0-2. 1M at m* N /rriN = 0.55 to about 1.39-1.63M 
at m* N /m N = 0.8. The mass splitting at m* N /m N « 0.65 for different K 
is about O.2M , increases to nearly O.35M at m* N /rriN ~ 0.75, and then 
again decreases to about O.25M at m* N /mN ~ 0.8. For comparison we 
also mark in Fig. [7] the maximum masses and the effective masses of sev- 
eral other RMF sets fitted to properties of nuclei, like TM1, NL3 or NL-SH 



271 142[ l25l . |43|, |44|. For the "model aupcp" , the maximum masses range from 



2.4 M at m* N /m N « 0.55 to 1.6 M at m* N /m N « 0.8. We also note that the 
results do not change substantially from the results obtained using "model 
aojp" when we include the a* meson (in "model aupa*4>"). Contrary to the 



statement made in [23j, |45[, that K as well as m* N "strongly influence the 
high-density behaviour of the equation of state of symmetric nuclear mat- 
ter" this being an "effect [which] shows up in neutron rich and neutron star 
matter", we find that the compression modulus has nearly no influence on 
the high-density behaviour of the nucleonic EoS whereas the effective mass 



has a stark impact. This finding is well known for the RMF model [191 . 146 
The effect of the compressibility on hyperonic stars is however remarkable, 
although this dep endence also varies considerably with the effective nucleon 



aep 

mass. In [23|, |45[, the investigated range of the effective nucleon mass was 



0.7-0.8 for the hyperonic case, in which K has a greater influence on the max- 
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imum mass than m* N /mN, which we also confirm in our study. However, we 
also extend our investigation for a broader range of m* N /m N , and find that 
this conclusion no longer holds for low effective nucleon masses. Also, at 
large m* N /m^ where K is influential the maximum masses of hyperonic stars 
do not comply with the limiting pulsar mass of 1.97±O.O4M . In "model 
acop" with hyperons this mass is only reached for low values of the effective 
nucleon mass m^/mjy < 0.6. Within the "model aujp(p" the mass constraint 
can only be achieved for m* N /m N < 0.7. In the case of nucleonic stars, for 
^jv/^Af ~ 0.78 the pulsar mass constraint cannot be reached if the com- 
pression modulus is below 225 MeV. 

To explore the whole range of K continuously, we plot lines of constant 
maximum mass of nucleonic stars in the m* N — K plane in Fig. [HJ The re- 
sults support our previous observation: the lines of constant maximum mass 
are nearly parallel to the x-axis, which means that the compression modulus 
does not influence the stiffness of the EoS at neutron star densities. Only 
for effective masses above m* N jmjq > 0.7 do the lines become slightly slanted 
and show that a higher compression modulus then slightly stiffens the EoS. 
We note, that in the investigated case of purely nucleonic neutron stars the 
observed mass of PSR J16 14-2230 requires an effective nucleon mass at sat- 
uration density of maximally m* N /m N < 0.82 at K — 300 MeV down to 
m* N /m N < 0.78 at K = 200 MeV. 

5. Summary 

The possibility of existence of hyperons in massive neutron stars has been 
investigated in this paper. Within hadronic models, we varied the coupling 
strengths of the scalar a meson to the E and 5 hyperons by changing the po- 
tential depths of these hyperons in nuclear matter. Meanwhile, we assumed 
SU(6) relations for the vector couplings and kept the A potential fixed. We 
found that for deeper potential depths the EoS becomes softer and the max- 
imum masses correspondingly smaller. Which of the two potentials has a 
larger impact on the stiffness of the EoS depends mainly on their relative 
values: if ui N) < U^ N) the Es will appear earlier (i.e. at lower densities) 
than the Es and therefore have the greater influence on the EoS. However, 
this cannot exactly be mirrored, since for ui N) > U^ N) the E~ might appear 
earliest of all Es and 5s, but as soon as the E~ appears, the number density 
of the E~ overtakes that of the E~ by far. The overall effect of varying the 
two potential depths is rather small (Am ^0.2M q over the probed range). 



19 




Figure 8: Lines of constant maximum mass of purely nuclear matter stars in the m* N — K 
plane. 
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Since the impact of the potential depths on the EoS is not enough to 
raise the maximum masses of neutron stars for "model aup" for GM1 above 
the limit of 1.97±O.O4M , we decided to increase the repulsive interaction 
amongst hyperons by including the vector <fr meson in the RMF model. The 
additional repulsion between hyperons results in a sufficiently stiff EoS ac- 
companied by an increase in the maximum masses of about Am ~O.2M . 
For the GM1 parameter set, the investigated area of the U^ N) - ui N) plane 
complies now to almost full extent with the mass of PSR J1614-2230. 

In the case of purely nucleonic stars we varied the effective mass of the nu- 
cleon as well as the compression modulus at nuclear saturation density which 
are input parameters for the RMF model. We found that the compression 
modulus has very little influence on the maximum mass of pure nucleonic 
neutron stars: in the range K = (200 — 300) MeV its influence is negligible 
for low effective masses m* N /rriN < 0.6, while its most prominent effect is 
for large effective masses m* N /rriN > 0.7 where it causes a mass difference of 
just A M < O.1M . The effective mass of the nucleon at saturation proved 
to be a very good parameter for obtaining large maximum masses. Over the 
investigated range m* N /m^ = 0.55 — 0.8 the maximum masses change about 
Am ~ 1.1M . The most massive stars (up to ~2.9M ) are obtained for low 
effective nucleon masses independently of the compression modulus, while for 
Tn* N /m,N > 0.78 we cannot reach the mass limit of 1.97±0.04M Q if the com- 
pression modulus is too low. For the hyperonic case, maximum masses varied 
from 2.0-2. 1M at m* N /m N = 0.55 to about 1.39-1.63M at m* N /m N = 0.8 
in the "model aup" , whereas in the "model aupcj)" the variation is from 2.4 
to 1.6 M & respectively. Hyperonic neutron stars in the "model aup" are only 
compatible with the new pulsar mass constraint for low values of the effective 
nucleon mass m* N /rriN < 0.6 MeV, while within the "model aup<p" the mass 
constraint can only be achieved for m* N /m^ < 0.7 MeV. We conclude that 
the nuclear compression modulus at saturation might be a good indicator for 
the stiffness of the nuclear matter EoS at low densities (around saturation), 
but the maximum mass is much more sensitive to the effective nucleon mass 
at saturation. This result is in accordance with the fact that the high density 
behaviour of the EoS in the RMF model is solely controlled by the effective 
nucleon mass at saturation density by virtue of the Hugenholtz-van Hove 
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As we have seen, a 2M star like PSR J1614-2230 helps to constrain 
the nuclear matter EoS. Then a possible 2.4M star - as has been recently 
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suggested to exist |47j - would pose even tighter constraints on the neutron 
star composition. Should this measurement be confirmed, we could rule out 
many hadronic models (e.g. TM1 or GM3) even in the limit of purely nu- 
cleonic stars. Furthermore, in most such models the appearance of hyperons 
together with the use of SU(6) relations for the vector coupling constants 



would then become impossible. In an associated paper [lif . we further in- 
vestigate the possibility of existence of hyperons in massive neutron stars, 
by questioning the assumption of SU(6) symmetry for the determination of 
vector meson-hyperon couplings. 
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